QCD sum rule calculation for the charmonium-like structures in the Jlip<p and J/if/oj invariant mass 
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Using the QCD sum rules we test if the charmonium-like structure 7(4274), observed in the J/t//<p invariant 
mass spectrum, can be described with a D J D t0 (2317) + h.c. molecular current with J pc = (T + . We consider the 
contributions of condensates up to dimension ten and we work at leading order in a s . We keep terms which are 
linear in the strange quark mass m s . The mass obtained for such state is m DsDs0 = (4.78 ± 0.54) GeV. We also 
consider a molecular 0~ + DD (2400) + h.c. current and we obtain m DDo = (4.55 ± 0.49) GeV. Our study shows 
that the newly observed 7(4274) in the J/if/(f> invariant mass spectrum can be, considering the uncertainties, 
described using a molecular charmonium current. 

PACS numbers: 14.40.Rt, 14.40.Lb, 11.55.Hx 



In the recent years, many new charmonium states were ob- 
served by BaBar, Belle and CDF Collaborations. There is 
growing evidence that at least some of these new states are 
non conventional cc states. In some cases the masses of these 
states are very close to the meson-meson threshold, like the 
X(3872) fl] and the Z + (4430) H. Therefore, a molecular 
interpretation for these states seems natural. Other possible 
interpretations for these states are tetraquarks, hybrid mesons, 
or threshold effects. 

Very recently the CDF Collaboration reported a further 
study of the structures in the J/iff(/> invariant mass, produced 
in exclusive B + — > J/if/<pK + decays. Besides confirming the 
7(4140) state fl] with a significance greater than 5cr, CDF 
also find evidence for a second structure with approximately 
3.1o" significance. The reported mass and width of this struc- 
ture are M = 4274.4+^(stat) MeV and V = 32.3+g|(stat) 
MeV JH. This new structure, refered as 7(4274) in ref. H], 
was interpreted as the S-wave D J D J o(2317) + h.c. molecular 
state. The authors of ref. JH] have also predicted a S-wave 
DD()(2400) +h.c. molecular state with a mass around 4.2 GeV, 
which they call as the cousin of 7(4274). This state is compat- 
ible with the enhancement structure around 4.2 GeV observed 
in the J/i//u> invariant mass spectrum from B decay |@]. 

These two pseudoscalar molecular states could be the ana- 
logue of the 7(4140) and 7(3930) (reported by the Belle Col- 
laboration |0] and confirmed bythe BaBar Collaboration 10]), 
that were interpreted, in ref. Q8J1 , as D* S D* S and D*D* scalar 
molecular states respectively. Some interpretations for the 
7(4140) can be found in refs. HUH]. 

Here we use the QCD sum rules (QCDSR) |j23i-|26ll, to 
check the suggestion made by the authors of ref. H . There- 
fore, we study the two-point function based on a D s D s q molec- 



ular current with J PC = 0~ + , to see if the new observed struc- 
ture, the 7(4274), can be interpreted as such molecular state. 
We also investigate the DDq molecular current. Previous cal- 
culations for the new charmonium states interpreted as molec- 
ular or tetraquark states can be found at lH3l l27l - l36ll . 

A possible D s D s q molecular current with J PC = 0~ + is 
given by 

j = —p [(Sa75Ca)(CbSb) + {c a J 5 S a )(s b C b )\ , (1) 

V2 

where a and b are color indices. 

The QCDSR approach is based on the two-point correlation 
function 



n(?) 



rf 4 * e ^<o|r[j(*)/(0)]|o>. 



(2) 



The sum rule is obtained by evaluating the correlation func- 
tion in Eq. (O in two ways: in the OPE side and in the phe- 
nomenological side. In the OPE side we work at leading or- 
der in a s in the operators, we consider the contributions from 
condensates up to dimension ten and we keep terms which are 
linear in the strange quark mass m s . In the phenomenological 
side, the correlation function is calculated by inserting inter- 
mediate states for the D s D s o molecular state. The coupling 
of the molecular state, M, to the current, j, in Eq. (HJ can be 
parametrized in terms of the parameter A 



(Q\j\M) = A. 



(3) 
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Although there is no one to one correspondence between 
the current and the state, since the current in Eq. (fl3 can be 
rewritten in terms of a sum over tetraquark type currents, by 
the use of the Fierz transformation, the parameter A, appear- 
ing in Eq. (O, gives a measure of the strength of the cou- 
pling between the current and the state. Besides, as shown in 
ref. 026I1 . in the Fierz transformation of a molecular current, 
each tetraquark component contributes with suppression fac- 
tors that originate from picking up the correct Dirac and color 
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indices. This means that if the physical state is a molecular 
state, it would be best to choose a molecular type of current 
so that it has a large overlap with the physical state. There- 
fore, if the sum rule gives a mass and width consistent with 
the physical values, we can infer that the physical state has a 
structure well represented by the chosen current. 

Using Eq. (0, the phenomenological side of Eq. (01 can be 
written as 



n phe "(q 2 ) 



q- 



f 

Jo 



ds'- 



s — q'- 



where the second term in the RHS of Eq.© denotes the con- 
tribution of the continuum of the states with the same quantum 
numbers as the current. As usual in the QCDSR method, it is 
assumed that the continuum contribution to the spectral den- 
sity, p ccn, (s) in Eq. ©, vanishes below a certain continuum 
threshold so. Above this threshold, it is given by the result 
obtained in the OPE side. Therefore, one uses the ansatz 031 



p co "'(s)=p OPE (s)®(s-s ). 



(5) 



In the OPE side the correlation function can be written as a 
dispersion relation: 



ds 



p UPE (s) 



s - q 



2 ' 



(6) 



where p 0PE (s) is given by the imaginary part of the correla- 
tion function: np 0PE (s) = Im[n 0/,£ (s)]. After transferring the 
continuum contribution to the OPE side, and after performing 
a Borel transform, the sum rule for the state described by a 
D s D s o pseudoscalar molecular current can be written as: 



,2 -m 
A e 



W M2 = ds e- s/M2 p 0PE {s) , (7) 



where 
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p OPE (s) = J]pW(s) 

D=0 



(8) 



with p [D1 representing the dimension-D condensates. To ex- 
tract the mass of the state we take the derivative of Eq. (|7} 
with respect to 1 /M 2 , and divide the result by Eq. (|7}: 



_ J4m, 
l D s D s0 ~ —FZE 



£\dse-^ 2 sp 0PE (s) 
£, dse-^-pOPE (s) ■ 



(9) 



The contributions to p 0PE (s), up to dimension-ten conden- 
sates, using factorization hypothesis, are given by: 
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where the integration limits are given by a ml „ = 
(1 - VI - 4/n?/s)/2, a max = (1 + -4m 2 c /s)/2, /?„„■„ = 
am 2 c l{sa — m 2 ), and we have used (sgcr.Gs) = m ( 2 (ss). For 
consistency, we have included the small contribution of the 
dimension-six condensate (g 3 G 3 ). We have also included 
the dimension-8 and dimension- 10 condensate contributions, 
related with the mixed condensate-quark condensate, gluon 
condensate squared, mixed condensate squared, four-quark 
condensate-gluon condensate and, three-gluon condensate- 
gluon condensate. 

For a consistent comparison with the results obtained for 
the other molecular states using the QCDSR approach, we 
have considered here the same values used for the quark 
masses and condensates as in refs. l27l - [36l l38tl : m c (m c ) = 
(1.23 + 0.05) GeV,m s = (0.13 + 0.03) GeV, (qq) = -(0.23 + 
0.03) 3 GeV 3 , <ss> = 0.8(qq), m\ = 0.8 GeV 2 , (g 2 G 2 ) = 
0.88 GeV 4 . For the three-gluon condensate we use (g 3 G i ) = 



0.045 GeV° 123 
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FIG. 1: The OPE convergence for the J pc = + , D s D sa molecule, 
in the region 2.2 < M 2 < 4.0 GeV 2 for yfa = 5.2 GeV. We plot 
the relative contributions starting with the perturbative contribution 
(solid line), and each other line represents the relative contribution 
after adding of one extra condensate in the expansion: + D = 3 
(dashed line), + D = 4 (dotted line), + D = 5 (solid line with circles), 
+ D = 6 (dashed line with squares), + D = 8 (dotted line with 
triangles), + D = 10 (solid line with triangles). 

The continuum threshold is a physical parameter that 
should be determined from the spectrum of the mesons. The 
value of the continuum threshold in the QCDSR approach is, 
in general, given as the value of the mass of the first excitated 
state squared. In some known cases, like the p and J/i//, the 
first excitated state has a mass approximately 0.5 GeV above 
the ground state mass. In the cases that one does not know 
the spectrum, one expects the continuum threshold to be ap- 
proximately the square of the mass of the state plus 0.5 GeV: 
so = (mo + 0.5 GeV) 2 . Therefore, to fix the continuum thresh- 
old range we extract the mass from the sum rule, for a given 
so, and accept such value of so if the obtained mass is in the 



range 0.4 GeV to 0.6 GeV smaller than yfso. Using this crite- 
rion, we obtain so in the range 5.1 < yfso < 5.3 GeV. 

The Borel window is determined by analysing the OPE con- 
vergence, the Borel stability and the pole contribution. To de- 
termine the minimum value of the Borel mass we impose that 
the contribution of the higher dimension condensate should be 
smaller than 10% of the total contribution: M 2 njn is such that 



OPE summed up dim n-1 (M 2 . ) 



total contribution (M~ h ) 



0.9. 



(11) 



In Fig. 03 we show the relative contribution of all the 
terms in the OPE side of the sum rule, in the region 2.2 < 
M 2 < 4.0 GeV 2 , for V^o = 5.2 GeV. From this figure 
we see that the contribution of the dimension- 10 condensate 
is smaller than 10% of the total contribution for values of 
M 2 > 2.4 GeV 2 , and that we have an excellent OPE con- 
vergence for M 2 > 2.4 GeV 2 . To have an idea of the impor- 
tance of the different terms in the OPE, we show, in Fig. [2] 
the contribution of each condensate. As we can see, the con- 
densates of dimension higher than six are, at least, one order 
of magnetude smaller than the perturbative contribution, in all 
considered Borel region. 
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FIG. 2: The OPE convergence for the J PC 



o- 



D s D s q molecule, 

in the region 2.2 < M 1 < 4.0 GeV z for y*o" = 5.2 GeV. We plot 
the contributions of all individual condensates in the OPE: the per- 
turbative contribution (solid line), {ss) contribution (dashed line), 
(g 2 G 2 ) contribution (dotted line), m^(ss) (solid line with cicles), 
(is) 2 (dashed line with circles), (g 3 G 3 ) (dotted line with circles), 
mj^ss) 2 (solid line with squares), {g 2 G 2 ) 2 (dashed line with squares), 
{m^(ss)) 2 (dotted line with squares), {ss) 2 (g 2 G 2 ) (dashed line with 
triangles), (g 2 G 2 )(g 3 G 3 ) (solid line with triangles). 

As commented above, the OPE convergence is very good 
in the Borel range 2.4 < M 2 < 4.0 GeV 2 . However, the 
Borel stability for the mass of the state is only good for 
M 2 > 2.8 GeV 2 , as can be seen through the solid line in 
Fig. [3] Therefore, we fix the lower value of M 2 in the sum 
rule window as M 2 . - 2.8 GeV 2 . 
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FIG. 3: The pseudoscalar meson mass, described with a D s D s o 
molecular current, as a function of the sum rule parameter (M 2 ) 
for s/sq = 5.2 GeV. The solid line shows the result obtained 
considering all contributions up to dimension- 10. The dashed 
and dotted lines show, respectively, the results obtained neglecting 
the contributions of the dimension-8 (Xg 2 G 2 ) 2 ) and dimension- 10 
«g 2 G 2 )<g 3 G 3 ), (g 2 G 2 )(ss) 2 ) gluon condensates. 

In Fig. [3] we also show, through the dashed and dotted 
lines, the result obtained if we neglect the contribution of 
dimension-8 and dimension-10 gluon condensates. We see 
that the contribution of the dimension-8 and -10 gluon con- 
densates ((g 2 G 2 ) 2 , (g 2 G 2 ){g 3 G 3 ), (g 2 G 2 ){ss) 2 ) are only im- 
portant in the region M 2 < 2.8 GeV 2 , which is not in our Borel 
window, due to the mass stability. Therefore, at least in this 
case, the contribution of higher dimension gluon condensates 
could be safely neglected. Besides, as can be seen in Fig. HI 
where we show the results for the mass for different values of 
sfso considering all condensate contributions up to dimension- 
10, the dependence of the mass on the OPE convergence is 
smaller than the its dependence on the continuum threshold 
parameter. 

To be able to extract, from the sum rule, information about 
the low-lying resonance, the pole contribution to the sum rule 
should be bigger than, or at least equal to, the continuum con- 
tribution. Since the continuum contribution increases with 
M 2 , due to the dominance of the perturbative contribution, 
we fix the maximum value of the Borel mass to be the one for 
which the pole contribution is equal to the continuum contri- 
bution. 

From Fig. [5] we see that for yso = 5.2 GeV, the pole 
contribution is bigger than the continuum contribution for 
M 2 < 3.66 GeV 2 . We show in Table I the values of M max 
for other values of yso". Although for yso =5.0 GeV there 
is still a small allowed Borel window, the difference between 
the obtained mass and the continuum threshold is very small 
(smaller than 0.2 GeV). Therefore, we do not consider values 
of v^o < 5.1 GeV. 

Table I: Upper limits in the Borel window for the 0~ + , 0,-0,0 



FIG. 4: The pseudoscalar meson mass, described with a D s D s q 
molecular current, as a function of the sum rule parameter (M 2 ) 
for yfso = 5.1 GeV (dotted line), yjjj = 5.2 GeV (solid line) and 
yfso = 5.3 GeV (dot-dashed line). 
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FIG. 5: The solid line shows the relative pole contribution (the pole 
contribution divided by the total, pole plus continuum, contribution) 
and the dashed line shows the relative continuum contribution for 
= 5.2 GeV. 

current obtained from the sum rule for different values of -\fso- 



V^(GeV) 


M,i. t (GeV 2 ) 


5.1 


3.43 


5.2 


3.66 


5.3 


3.90 



To estimate the dependence of our results with the values 
of the quark masses and condensates, we fix yfso = 5.2 GeV 
and vary the other parameters in the ranges: m c = (1.23 ± 
0.05) GeV, m s = (0.13 + 0.03) GeV, (qq) = -(0.23 ± 
0.03) 3 GeV 3 , ml = (0.8 + 0.1) GeV 2 . In our calculation we 
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have assumed the factorization hypothesis. However, it is im- 
portant to check how a violation of the factorization hypoth- 
esis would modify our results. To do that we multiply the 
contribution of the four-quark condensates of D = 6, 8 and 
D — 10 in Eq. (TTOb by a factor K, and we vary K in the range 
0.5 < K < 2. The dependence of our results with all the 
variations mentioned above is show in Table II. 

Table II: Values obtained for m DsDs0 , in the Borel window 
3.0 < M 2 < 3.65 GeV 2 , when the parameters vary in the ranges 
showed. 



parameter 


'"D,,D,, (GeV) 


m c = (1.23 ±0.05) GeV 


4.76 ± 0.07 


m s = (0.13 ±0.03) GeV 


4.76 ± 0.06 


<^> = -(0.23 ± 0.03) 3 GeV J 


4.89 ±0.27 


m 2 = (0.8 ±0.1) GeV 2 


4.753 ± 0.003 


0.5 < K < 2 


4.80 ±0.11 



Taking into account the uncertainties given above and the 
uncertainties due to the continuum threshold parameter and 
due to the OPE convergence, we finally arrive at 



m DsDM = (4.78 + 0.54) GeV, 



(12) 



which, considering the error, is still in agreement with the 
mass of the newly observed structure F(4274). 

One can also deduce, from Eq. (0, the parameter A defined 
in Eq. 0]). We get: 



(6.0 + 3.9) x 10~ 2 GeV 5 . 



(13) 



This number is of the same order as the current-state coupling 
obtained in ref. lfl3ll . where the J PC = ++ D*D* molecular 
current was considered to describe the 7(4140): 

A D . D . = (4.22 + 0.83) x 10~ 2 GeV 5 . (14) 

Therefore, we can conclude that the state can be well repre- 
sented by the J PC = 0~ + D s D s q molecular current 

To obtain results for the DDq molecular current with J PC = 
0~ + , we only have to take m s = and (ss) = (qq) in Eq. ( flOl . 
As can be seen by Fig. [6] the OPE convergence in this case 
is also very good for M 2 > 2.4 GeV 2 . Therefore to fix the 
minimum value of the Borel parameter, we will consider the 
Borel stability of the obtained mass. For this we show, in 
Fig.|7]the results for the mass of the state described by a DDo 
pseudoscalar molecular current, for different values of yfsQ. 
We see that for M 2 > 2.7 GeV 2 we get a good Borel stability. 
Therefore we fix M 2 min = 2.7 GeV 2 . 

In Table III we give the values of M max for the considered 
values of yl7j. 

Table III: Upper limits in the Borel window for the 0~ + , DD 
current obtained from the sum rule for different values of -Jsq. 



V^(GeV) 


Ml. t (GeV 2 ) 


4.9 


3.19 


5.0 


3.40 


5.1 


3.61 



Taking into account the uncertainties due to the quark 
masses, condensates, continuum threshold parameter and 
OPE convergence, we finally arrive at 
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FIG. 6: Same as Fig. 1 for the DDo current for -\fso = 5.0 GeV: per- 
turbative contribution (long-dashed line), relative contribution after 
adding of one extra condensate in the expansion: + D = 4 (dotted 
line), + D = 6 (dashed line with squares), + D = 8 (dotted line with 
triangles), + D = 10 (solid line with triangles). 
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FIG. 7: Same as Fig.f4]for the DDo molecular current, for yfso = 4.9 
GeV (dotted line), y^o = 5.0 GeV (solid line) and y^o = 5.1 GeV 
(dot-dashed line). 



which, although a little bigger than the prediction in ref. JH 
for a S-wave DDq molecular state, is still in agreement with it, 
considering the error. It is interesting to notive that the result 
in Eq. ( fT3T > is in a excellent agreement with the result obtained 
in ref. [40], where different tetraquark currents were used to 
study J PC = and 0~ + charmonium-like states. For the 
parameter A we get: 



*dd = (5.4 + 3.9) x 10 ~ GeV 5 



(16) 



m DDo = (4.55 ± 0.49) GeV, 



(15) 



The mass we have obtained for the DDq molecular state 
is approximately two hundred MeV below than the value ob- 
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tained for the similar strange state. This is very different from 
the results obtained in ref. IU3I1 where the J PC = ++ D*D* S and 
D*D* molecular currents were considered. In the case of the 
scalar molecular currents, the difference between the masses 
of the strange and non-strange states was consistent with zero. 

In conclusion, the newly observed structure F(4274) in the 
J/ijjtp invariant mass spectrum can be, considering the errors, 
interpreted as the S-wave D s D S Q+h.c. molecular charmonium, 
in agreement with the findings in ref. JH, where a dynami- 
cal study of the system, composed of the pseudoscalar and 
scalar charmed mesons, was done. In the case of the S-wave 
DD() + h.c. molecular current, which was called as the cousin 
of 7(4274) in ref. 0], the QCDSR results are consistent with 
the enhancement structure around 4.2 GeV in the J/if/u> invari- 



ant mass spectrum from B decay JfJ. 
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